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Abstract 

Using e + e~ — > D*~D^~ data collected near the peak D s production energy, E cm = 4170 MeV, 
with the CLEO-c detector, we present the first observation of the decay Df — > pn. We measure 
a branching fraction B(Df — > pn) = (1.30 ± 0.36to'ie) x 10 3 - This is the first observation of a 
charmed meson decaying into a baryon-antibaryon final state. 
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Of the three ground state charmed mesons, only the Df is massive enough to decay 
to a baryon-antibaryon pair. Even before the discovery of the D+, a search for the decay 
Df — ► pn was suggested |l| as a "smoking gun" for decays proceeding via annihilation 
through a virtual W + , and a prediction was made that the branching fraction would be 
~ 1% if the annihilation mechanism dominated decays. In the intervening period it 
has become clear that the annihilation diagram contributes to, but does not dominate, 
decays, and has been studied in purely leptonic decays such as Df —> fiv Q and — > tv 
(3|. However, although the theoretical study of Df — > pn is complicated by final state 
interactions, it still has a unique role to play in the understanding of charmed meson decays. 

Finding decay modes that include an (anti-)neutron is particularly challenging. CLEO-c 
is the first detector to have a large dataset (~ 325 pb _1 ) of e + e~ annihilation events taken 
at a center-of-mass energy of around 4170 MeV. At this energy, there is a substantial cross- 
section (»s 1 nb) for the reaction e + e~ — ► D*~Df Q]. Using the knowledge of center-of-mass 
energy and momentum, we use a missing mass method to find the (anti-) neutron and so 
do not depend upon its interaction in the detector. Therefore, in this paper, mention of a 
particular decay mode implies the use of the charge conjugate decay mode also. 

The CLEO-c detector [5j is designed to measure the momenta of charged particles which 
curve in a 1.0 T solenoidal magnetic field, and identify them using specific energy loss 
(dE/dx) and Cherenkov imaging (RICH). Photons are detected, and their energy measured, 
using a Csl calorimeter. Our analysis procedure has much in common with that used in the 
measurement of B(D+ — > fiu) [2J. Here, we fully reconstruct one D~ as a "tag," reconstruct 
a transition photon from a D* decay, and identify and measure the momentum of a proton. 
We can then reconstruct the missing mass of the event and look for a peak at the anti-neutron 
mass. 

The D~ tags are found in the eight modes: K + K~n~ , K®K~ , rjir", r)'ir~, 4>p~ , n~ii + i{~ , 
K*~K*° and i]p~. Track selection, particle identification and definition of resonances are 
similar to those in our previous publication 0], with one important exception; each of the 
charged tracks in the D~ tag is required to have a dE/dx measurement more than three 
standard deviations, a, away from that expected for a proton. All the D~ tags are required 
to have momentum consistent with coming from the two-body production D*D S . Figure 1 
shows the measured mass of the tag candidate minus the nominal D~ mass, divided by 
the resolution of the candidate tag's mode. The fit shown is a unit Gaussian centered at 
zero. Those candidates within 2.5 of the peak are kept as D~ tags. According to the fit they 
comprise 27700 real D~ mesons and 64900 background combinations. Those combinations in 
the regions ±(3.5 — 6.0) in this plot will be used as a check on the combinatorial background 
to our final signal. We kinematically constrain the mass of the tags in the signal region to 
the known mass of the D~ . To be consistent, each sideband tag is kinematically constrained 
to the center of its sideband. 

We now add a 7 that satisfies our shower shape requirements to the D~ tag. Using the 
four- momenta of the 7 candidate, ^(7), and the D~ tag, we calculate the four-momentum 
of the .D+ using the equation p(Df) = Pbeam — p(D~) — p(j). The four-momentum of 
the beam, Pbeam takes into account the small crossing angle of the CESR beams. The 
missing mass squared (MM 2 ) distribution for good tag events is shown in Fig. 2. It shows 
a peak at M 2 (Df) corresponding to D*D S production. This is fit to a signal shape of a 
Crystal Ball function 0] with fixed tail parameters derived from Monte Carlo simulation, 
together with a fifth order polynomial background function. We select those events with 
MM 2 values between 3.779 and 3.976 GeV 2 . This is a loose requirement, with most of the loss 
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FIG. 1: The reconstructed mass minus the known D s mass, divided by the detector resolution, 
for all eight modes of D s tags reconstructed. The fit shown is a unit Gaussian centered at zero, 
together with a second order polynomial background function. 

in efficiency due to initial state radiation of the beam, which smears the MM 2 to artificially 
high values. According to the fit, the yield of D~ r yD^ candidates in this range is 16955 
above a background of 63170. This yield will be the denominator in our final branching 
fraction calculation. 

We next select our proton candidate. Monte Carlo simulation shows that all protons 
from this decay mode will have momenta in the range 150-550 MeV/c. This is below the 
momentum range for the RICH detector to identify protons, but well suited to identification 
by dE/dx. We require that this measurement be within 3a of that expected for a proton, 
and greater than 3cr from that expected for a kaon or a pion. The overall proton efficiency 
is determined by Monte Carlo simulation to be 75%, with the efficiency lowest at the lower 
proton momenta. 

We may now calculate the missing 4-momentum in the event, equal to the expression 
Pbeam—p(D~) —p(j) —pproton, and can thus calculate the missing mass of the event. However, 
we have further kinematic constraints we can impose which allow us to improve the missing 
mass resolution and reject combinatorial background. We do not know a priori if the photon 
is due to the transition D*~ — > D~(tag)j, or D* + — > D+ (signal)^ . We perform kinematic 
fits with each assumption, and choose between the two based on the x 2 values of the two 
fits. First, we add the photon to the D~ tag to form a D*~ candidate, and constrain 
the momentum of this D*~ candidate to that calculated from the two body production 
e + e~ — >• D*~D+. We then constrain the mass difference M(D*~) — M(Dj) to its nominal 
value. Alternatively, we constrain the Dj tag itself to the momentum calculated assuming 
the two-body production e + e~ — >■ D S ~D* + , then combine the proton with the missing mass 
of the event to make a D s signal candidate, add the photon, and constrain the M(D* + — Df) 
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FIG. 2: The missing mass squared from events with a reconstructed 7 and D s {tag). The fit is to 
a signal shape of a Crystal Ball function [6] with fixed tail parameters, together with a fifth order 
polynomial background function. 



mass difference. We choose the scheme with the lowest total x 2 value; in Monte Carlo 
simulation we find that we assign the photon to the correct D s greater than 95% of the time. 
The kinematic constraints on the detected particles improve the resolution in missing mass 
by around a factor of two, whichever D s the photon is combined with. Furthermore, we can 
place cuts on the x 2 of the kinematic constraints to reject combinatorial background. In the 
case of the momentum constraint we require x 2 < 9> an d in the case of the mass-difference 
constraint we require x 2 < 4; with each constraint there is one degree of freedom. The 
requirement is looser for the momentum constraint because initial state radiation produces 
a tail in the momentum distribution. 

The transition photon in the event has an energy in the laboratory of 110-180 MeV. In this 
energy range there is the possibility of background clusters passing all the requirements for 
being a photon. Such background photons are particularly prevalent in events which contain 
anti-baryons as they frequently interact with the detector and give "split-off" clusters, often 
far from the impact point of the particle in the Csl calorimeter. Occasionally an event may 
survive all the above requirements while having more than one photon candidate. If so, 
we select the photon candidate that produces the lowest combined x 2 m the kinematic fit. 
Background photons also influence the signal shape which we determine using Monte Carlo 
simulation. This shape is well described by a core Gaussian function of a ~ 4 MeV centered 
at the neutron mass, together with a second, offset, Gaussian of width o « 38 MeV and 
containing « 12% of the signal. This second Gaussian is due to events where we have used 
an incorrect photon candidate. 

Figure 3 shows the missing mass distribution for the events after all requirements and 
kinematic fitting, and contains thirteen events. These are the only events in the missing- 
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FIG. 3: The missing mass in the event after all requirements and kinematic fitting has been 
performed. The fit is described in the text. 



mass range 600-1100 MeV. The plot is well fit using a likelihood fit to the signal shape 
described above, so we take our signal yield to be the 13.0 ±3.6 events observed. Repeating 
the analysis using sidebands to the D~ as described above, gives three events in the missing- 
mass range 600-1100 MeV, none of which are in the signal region of 900-980 MeV. We divide 
this yield of 13 by the number of Df decays we have detected, and correct for the efficiencies 
of requirements placed on the fit \ 2 an d proton reconstruction and identification. This gives 
a branching fraction of (1.30 ±0.36) x 10~ 3 , where the error shown is the statistical error in 
the signal yield only. 

We have performed many checks to ensure that our analysis is not biased towards ob- 
taining events only in the signal region. To enter the final signal plot, an event must have 
a proton candidate of a momentum which happens to be well matched to the capabilities 
of the dE/dx system, and thus background from incorrectly identified proton candidates is 
negligible. This, in turn, means that background from other charm events is negligible. In 
order to check that the sidebands are a reasonable representation of combinatorial back- 
ground, we generated a large sample of uds continuum Monte Carlo events with final states 
of a proton, an anti-neutron, a photon, and the decay products of a D~. Very few of these 
events passed all the kinematic requirements, and among those that did, there were as many 
anti- neutron candidates in the sideband plot as the signal plot. 

As a check on the efficiency of the selection criteria imposed on the quality of the kine- 
matic fit, we repeated the analysis looking for the decay D£ — > K + K°, where we consider 
the K° to be the missing-mass analog of the antineutron. We reproduce very well the mea- 
surement of this branching fraction found by more direct means [7] , and from the comparison 
of the efficiencies of the kinematic constraint cuts in data and Monte Carlo, we find a sys- 
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tematic uncertainty in the efficiency of ±4%. The uncertainty in the number of Df, the 
denominator in our branching fraction calculation, is systematically limited, and by looking 
at the variation of the yield using a variety of different signal and background functions, 
estimated to be ±5%. One possible source of systematic uncertainty concerns the signal 
shape, the tail of which is dominated by split-offs from antibaryon interactions which are 
one of the hardest processes to reliably simulate. We conservatively assign a ±6% uncer- 
tainty to account for any mis-modelling of this process; this number corresponds to getting 
the number of events from these fake photons incorrect by ±50%. Lastly, we note that 
although there is no evidence of any background in the final plot, we cannot assume that 
it is strictly zero. As our best estimate of the background is events, which corresponds 
to an upper limit (1 a) of 1.1 events, we therefore introduce a -8.5% uncertainty in the 
final branching fractions. Combining these systematic uncertainties in quadrature, produces 
a total systematic uncertainty of i?2% in the branching fraction, much smaller than the 
statistical error. 

In conclusion, we report the first observation of the decay — > pn with a signal of 13 
events and a background consistent with zero. We measure the branching fraction B(D+ — > 
pn) = (1.30 ± 0.36io;ig) x 10~ 3 . This is the first observation of a charm meson decaying 
into baryon-antibaryon pair. The two-body decay observed here is the only one allowed 
kinematically. The actual decay process is suspected to be related to annihilation, which is 
also responsible for purely leptonic decays. Relating this baryonic decay rate to the leptonic 
rate should provide important clues as to how baryons are produced in hadronic interactions. 
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